Tunable thermal conductivity in defect engineered nanowires at low temperatures 
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We measure the thermal conductivity (k) of individual InAs nanowires (NWs) , and find that it is 
3 orders of magnitude smaller than the bulk value in the temperature range of 10 to 50 K. We argue 
that the low k arises from the strong localization of phonons in the random superlattice of twin- 
defects oriented perpendicular to the axis of the NW. We observe significant electronic contribution 
arising from the surface accumulation layer which gives rise to tunability of k with the application 
of electrostatic gate and magnetic field. Our devices and measurements of k at different carrier 
concentrations and magnetic field without introducing structural defects, offer a means to study 
new aspects of nanoscale thermal transport. 



Thermal transport measurements on semiconducting 
nanowires (NWs) have attracted a lot of attention in 
the last few years. Measurements of thermal trans- 
port in nanostructures are important as they provide 
a platform to test the existing descriptions of phonons 
in confined structures and across complex interfaces^, 
and have the potential to result in technological appli- 
cations as thermoelectric system^E2. Different materials 
are benchmarked using the thermoelectric figure of merit 
= 2-T where S is the Seebeck coefficient, p is the 
electrical resistivity, k is the thermal conductivity and 
T is the absolute temperature. ZT can be increased by 
appropriate engineering of nanostructures. One of the 
ways to increase ZT is by reducing k without degrading 
its electrical conductivity and Seebeck coefficients. As 
a result, an ideal thermoelectric material is a glass for 
phonons and ordered for electronic transport. Recently, 
several theoretical models as well as experimental stud- 
ies have been carried out in different semiconductor NWs 
like Si, Ge, Bi 2 Te 3 etcPH. It is found that for Si NWs 
the value of k is reduced by two orders of magnitud e! 5 * 8 * 9 ! 
compared to bulk values by tuning the roughness of the 
surface. III-V semiconductors are also known to be good 
thermoelectric materials, and theoretical studies suggest 
that InSb and InAs NWs are good candidates for bet- 
ter ZT^. InAs NWs have been studied extensively to 
probe their charge and spin-transport!^^. An aspect 
of InAs NWs that makes s tudyi ng their thermal trans- 
port, hitherto little explore d 14 * 15 * , of interest is the ability 
to tune the density of twin defects and polytypes along 
its length by varying growth parameter d 16 " 20 *. Exploit- 
ing this control over crystal structure can help synthe- 
size defect-engineered NWs, whose lattice has aperiodic 
array of twins along its length that modify phonon be- 
havior, without significantly compromising their electri- 
cal properties^. Such NWs satisfy the key criteria for 
a good thermoelectric material - localization of phonons 
without localizing electrons. 

In this work we explore the k of suspended InAs NW 
field effect transistors (FETs); the NWs have high density 




-15 -10 -5 5 10 15 
Gate Voltage (V) 



30 
T(K) 



FIG. 1. (color online) (a) Schematic of the four probe sus- 
pended NW device and the circuit used for measuring elec- 
trical signals for both electrical and thermal transport, (b) 
A tilted scanning electron microscope image of a device. The 
scale bar corresponds to 1 rim. (c) Resistance as a function 
of gate voltage V 9 at bath temperature of 30K for a typical 
device, (d) Resistance R as a function of temperature for a 
device at three different gate voltages V 9 . The solid curves 
show linear fits to the data. 



of aperiodic twins and polytypes perpendicular to their 
ax j^i6H18llo] ranc 5om nature of defects suppresses n p h 
(the phonon contribution to k) resulting in reduction of 
thermal conductivity by 3 orders of magnitude from bulk. 
Our method also has an advantage which is to study k 
as a function of carrier concentration without introduc- 
ing additional impurities. We observe significant change 
in the thermal conductivity with the application of gate 
voltage suggesting finite electronic contribution, and we 
propose that this contribution comes from the conduc- 
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tion electrons at the surface accumulation layer which is 
inherent to low band gap semiconductors like InAs and 
I n ]\J2iH23] -j> con fj rm the role of electronic contribution 
we apply magnetic field and see considerable changes in 
the thermal conductivity. 

The InAs NWs used to fabricate the devices in this 
work were grown in a metal organic chemical vapour de- 
position (MOCVD) system using the vapor-liquid-solid 
(VLS) technique on a < 111 > B oriented GaAs sub- 
strate. NWs of 100 - 200 nm diameter were used in 
this experiment. 10-12 /im long NWs were chosen for 
making four probe suspended devices on a 300 nm thick, 
SiC>2 coated degenerately doped Si wafer, which serves 
as the back gate (growth and device fabrication details 
are given in supplementary Sl^J). Fig. [ija) shows the 
schematic of the device and the circuit for electrical and 
thermal measurements. Fig. |TJt)) shows a scanning elec- 
tron microscope (SEM) image of a device. The gold elec- 
trodes provide mechanical support and Ohmic contacts 
with negligible contact resistance to the suspended NW 
besides thermally anchoring it at four points to the bath 
temperature Tj,. The electrical properties of the device 
are characterized as a function of gate voltage (V g ) and 
temperature before detailed thermal transport measure- 
ments. Fig.jljc) shows the resistance (R= ) of a typ- 
ical device as a function of V g at a fixed value of T&. 
The n-type field effect transistor (FET) behavior of the 
suspended NW devices is clearly seen. Fig. [ljd) shows 
the evolution of R as a function of temperature, T, and 
also the corresponding linear fit in the temperature range 
10-50 K. 

High resolution transmission electron microscopy 
(TEM) studies were used to characterize the microstruc- 
ture of our twin-engineered NWs. Fig.[2|a) shows a large 
scale TEM image of a NW where we observe an aperi- 
odic arrangement of the twin superlattice perpendicular 
to its axis. The fact that the lattice has defects but lacks 
a long-range periodicity is reflected in the Fourier ampli- 
tude (seen in Fig.[5]jb)) that has no sharp peaks along the 
length of the NW. Figs, [gfc) fc (d) show a detailed view 
of the twins and polytypej^ * 18 * 20 ! The typical spacing 
of these defects is ~3-5 nm. Once electrical and struc- 
tural characterization is done, we perform simultaneous 
thermal and electrical measurements. 

One method of measuring thermal properties of nanos- 
tructures uses a modification of the 3w technique^ for 
k and specific heat measurement. The modified 3w 
techniqueP^H re lies on fabricating suspended four probe 
devices with the conducting NW serving both as a heater 
and a thermometer. The magnitud e of V (3u>) can be re- 
lated to k (in the limit ujj — > ) a d 26 ' 27 1 
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FIG. 2. (color online) (a) TEM image of an InAs NW shows 
twin boundaries perpendicular to the axis of the NW. The 
scale bar is 50 nm. (b) Fourier amplitude of the intensity 
modulation along the length of the NW that shows no sharp 
peaks, (c) A close-up image showing the twin interfaces and 
the resulting zig-zag surface roughness of the NWs. The scale 
bar is 10 nm. (d) TEM image showing the atomic fringes and 
the orientation of the twin boundaries and polytypes. The 
scale bar is 2 nm. 



where u) is the angular frequency of the AC current, I(u>) 
is its amplitude, R and R' are the resistance and the 
derivative of resistance with respect to temperature re- 
spectively, L is the length of the NW between the two 



middle probes and A, its cross sect ion J2S12ZI. For our de- 
vices, L is 2 — 4/im and the thermal time constant 7 can 
be estimated to be ~ 10 _5 s (details of phase informa- 
tion of 3cj signal in the limit of W7 — > are given in the 
supplementary S32^). 

To measure k at each temperature value, we vary the 
current (I(w)) and obtain k from the fitting of Equation!!] 
to the measured V(3w) as a function of injected I(wJT 
Fig. [3^b) shows the evolution of V(3w) with I(w) along 
with the cubic power-law relationship that is essential to 
ensure that the 3lo technique is applicable within the ap- 
plied current regime (to eliminate other source that can 
possibly show cubic scaling see supplementary S^P^J). As 
we vary the temperature, we also record the resistance. 
We determine R' from the derivative to the linear fit to 
the R-T plot (see Fig.Qd)). The variation of R with T 
becomes small and non-monotonic above 50-60 K this is 
the reason we choose the temperature range of 10-50 K in 
our experiment. The phase of the 3uj signal is measured 
to be zero at all temperatures (see supplementary Fig. 
grj24]j xhis procedure of fitting the power law given by 
Equation [T] at each instance of parameter space allows 
us to extract k as a function of V ff and temperature. 
We have calibrated our technique using lithographically 
fabricated 170 nm wide 75 nm thick gold NWs in the 
temperature range of our measurement and obtained a 
Lorenz number around 2.4 x 10~ 8 WCl K~ 2 (see supple- 
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FIG. 3. (color online) (a) The temperature profile along the 
length L of the NW with the red shaded area indicating the 
amplitude of oscillations, (b) The measured V(3cj) signal as a 
function of the excitation current I(w) for three different tem- 
peratures T. The solid line shows the fit to the cubic power 
law (Equation [TJ . (c) Plot of k as a function of temperature 
for two different gate voltages, (d) Measured k as a function 
of temperature (dotted lines) together with electronic contri- 
bution re e ; (solid lines), (e) Measured k for the same device 
at three different magnetic fields as a function of the temper- 
ature at zero gate voltage. 

mentary S3^l) which calibrates the measurement tech- 
nique we have used. 

Fig. |3jc) shows the plot of measured function of 

temperature at two values of Y g . The overall magnitude 
of k is about 1000 times smaller than the bulk values^ 
(see supplementary Fig. S1<P^J). This is a significant 
reduction in k and the measured value is similar to the 
minimum thermal conductivity model of Cahill et a/P^ 
(CWP model) based on the assumption that the material 
is highly disordered. Similar reduction in k, beyond the 
CWP model, was also observed in disordered WSe 2 which 
is a layered structure^. Later we return to a model to 
understand this better. Similar evolution of k is also 
observed at different values of Y g in all the measured 
devices (see supplementary Fig. S 10-PP^). 

Now we consider the effect of Y g on k. It is clear that 
the electrical conductance of the semiconducting NW 
FET can be modified as a function of Y g (Fig. fife)). 
However, the effect of electron density variation is not 
expected to affect k in most semiconducting materials, 
as phonons are the dominant carriers of heatP. Di- 



rect calculation from electrical conductivity using the 
Wiedemann-Franz (WF) law gives a small value of the 
electronic contribution (/t e ;) which does not explain the 
variation we see with Y g (see supplementary S(P3). To 
understand the origin of large n e i we have to first un- 
derstand the material property of the InAs wire. It is 
well known that in InAs there is an accumulation of 
charge at the surface. The resultant downward band 
bending pins the surface Fermi level above the conduc- 
tion band minimuirPIHsII. In a geometry where the sur- 
face to volume ratio is large, such as NWs, this surface 
accumulation charge (SAC) plays a major role in the 
electrical transport. The SAC takes part in the elec- 
trical conduction and scatters due to roughness at sur- 
face, which is also the reason for low electron mobil- 
ity in InAs NWs compared to bulkJ2E2 We have pas- 
sivated our devices using (NH4)S X but this passivation 
leads to further enhancement to the surface charge den- 
sity (~ 10 12 cm _2 jp3 and consequently excellent Ohmic 
contacts. Fig. [ijc) shows that our device is always in 
conducting state throughout the gate range with low re- 
sistance value, strongly suggesting the presence of the 
accumulation layer. To estimate the n e i we have to con- 
sider the thickness of the layer of this SAC since it will 
reduce the effective cross sectional area through which 
heat is transported by electrons. The electrons in the 
bulk and in the accumulation layer can be thought of 
two parallel channels conducting heat and we can define 
a effective layer thickness (t) to find out the actual K e i 
from the WF law, 

^ = ^-V-^} {0<t - r) (2) 

where I and r are length and radius of the wire respec- 
tively, R is the electrical resistance at temperature Tb and 
L = 2.45 x 10~ 8 Wfl K~ 2 is Lorenz number. We can 
get a rough estimate of t from above equation by assum- 
ing that at lowest temperature the phonon contribution 
(n p h) is small and K e i ~ n. For the NW of length 2.6 
/jm and radius 70 nm we get t from our data as ^20 nm 
and ~10 nm for applied gate voltages of -10V and +10V 
respectively. Now we use Equation 2 to find n e i for all 
Tb, the result is plotted as solid lines in Fig. Ifd). This 
estimated K e i from a simplified picture captures the vari- 
ation we see in our data with V g . Another confirmation 
about this non negligible K e i comes from our magnetic 
field measurement. One way to suppres s the electronic 
contribution is to apply a magnetic field- 2 -^, as it will 
reduce the electronic mean free path for heat transport 
by deflecting the semiclassical trajectory of electrons. We 
have done measurement in magnetic field up to 6T and 
see considerable reduction of thermal conductivity (see in 
Fig.fJJe)). Magnetic field can affect phonons via electron- 
phonon interaction, the amount of change in tt p h however, 
is much smaller^ than the effect we see in our measure- 
ment. 

The total k is small and it is important to understand 
the origin of reduction in K p h- It is clear from the struc- 
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FIG. 4. (color online) The plot of measured function 
of temperature of two different NW samples at V 9 , along 
with, the predicted values from two models to understand 
the reduction in n p h for InAs NWs. Based on CWP model of 
Cahill et a/P, and HP model, of Hopkins et dP. 



tural analysis in Fig. [2] that our NWs have a high density 
of twins. The twin defects introduce additional inter- 
faces and also cause the surface to be rough, with typical 
roughness ~2-3 nm. The heat conduction in bulk is dom- 
inated by phonons whose mean free path \ m f P is much 
larger than the spacing between defect planes throughout 
the temperature range of our measurements (see supple- 
mentary Fig. SlfP^J). The average number of defect inter- 
faces is ~500 per /im - this makes the problem challeng- 
ing to model comprehensively. However, previous exper- 
iments on orderecPSEH, quasiperiodic superl attice f^ and 
calculations with modulation of superlattic d 38 * 39 ! on sim- 



ilar lengthscale << X m f P suggest that suppression of k 
takes place in such structures. Our devices have aperi- 
odic interfaces (see Fig. [2]Ja)) and they are likely to at- 
tenuate a wide band of phonons. Additionally, modeling 
indicates that heat transport across grain boundaries can 
offer significant thermal impedance^. In order to under- 
stand the magnitude of the observed k, we use a model by 
Hopkins et alW* (denoted by HP model) which explicitly 
takes into account the role of interfaces for calculating 
minimum k in layered materials. We have calculated the 
predicted values of the model of Cahill et al. ^ (de- 
noted by CWP model) and the HP model for InAs (see 
supplementary S£p^). There is a range of k measured for 
devices made from different NW samples (see supplemen- 
tary SEP2); data from two such devices at Y g is shown 
in Fig. [4j Without introducing any free parameters, the 
predictions of the CWP and HP models result in values 
that are comparable to the range in k measured in our 
experiments. 

In summary, we demonstrate that thermal conductiv- 
ity of InAs NWs with twin defects and polytypes is 1000 
times smaller than bulk values which will make our de- 
vices suitable candidates for good thermoelectric materi- 
als. We show the electronic contribution becoming signif- 
icant due to presence of the surface charge accumulation. 
The ability to tunc thermal conductivity electrostatically 
can lead to novel ways of probing heat transport in defect 
engineered nanostructures, and possible device applica- 
tions. 
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